GROWTH RATES AND EXPLOSIONS IN SANDPILES 
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Abstract. We study the abelian sandpile growth model, where n 
particles are added at the origin on a stable background configu- 
ration in Z^. Any site with at least 2d particles then topples by 
sending one particle to each neighbor. We find that with constant 
background height h < 2d — 2, the diameter of the set of sites 
that topple has order n^^'^. This was previously known only for 
h < d. Our proof uses a strong form of the least action principle 
for sandpiles, and a novel method of background modification. 

We can extend this diameter bound to certain backgrounds in 
which an arbitrarily high fraction of sites have height 2d ~ I. On 
the other hand, we show that if the background height 2rf — 2 is 
augmented by 1 at an arbitrarily small fraction of sites chosen in- 
dependently at random, then adding finitely many particles creates 
an explosion (a sandpile that never stabilizes). 



1. Introduction 

In this paper we consider the abehan sandpile model as a growth 
model in the integer lattice Z*^. The model starts from a stable back- 
ground configuration in which each site x has a pile of (j{x) < 2d — 1 
particles. To this background, n particles are added at the origin. 
Typically, n is large. We stabilize this configuration by toppling every 
unstable site; that is, every site with at least 2d particles gives one 
particle to each of its neighbors, until there are no more unstable sites. 
For more information on the abelian sandpile model, also known as the 
chip-firing game, see [1, 2, 3, 11]. 

To keep things simple in this introduction, we will enumerate the 
sites in (say, in order of increasing distance from the origin, breaking 
ties arbitrarily) and perform topplings one by one in discrete time: 
At each time step, if there are any unstable sites, then the smallest 
unstable site topples. All of our results hold also for the more general 
toppling procedures discussed in section 2. 
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Figure 1. Left: Stable sandpile of ri = 2 ■ 10^ particles 
in 1? on background height h = 2. Right: Sandpile of 
n = 15000 particles in on background height 3, except 
every fifth row and column has background height 2. In 
both cases, the set T„ is a square. Color scheme: sites 
colored blue have 3 particles, turquoise 2 particles, yel- 
low 1 particle, red particles. 



Let T„ = Tn^d,a be the set of sites that topple (Figure 1). Since these 
sets are nested, Ti C T2 C . . ., it is natural to view them as a growth 
model, with n playing the role of a time parameter. We distinguish 
between two extreme cases. If T„ is finite for all n, we say that a is 
robust. In this case we are interested in the growth rate, i.e. in how the 
diameter of T„ grows with n. 

At the other extreme, if T„ = Z*^ for some n, then every site topples 
infinitely often. Otherwise, some site x E Z'^ must finish toppling 
before all of its neighbors do; since each neighbor topples at least once 
after x finishes toppling, x receives 2d additional particles and must 
topple again. 

If Tn = Z'^ for some n, then we say that a is explosive, and it is 
exploding when the n particles are added. (In [7], the term 'not stabi- 
lizable' was used for 'exploding,' and 'metastable' for 'explosive.') The 
simplest example of an explosive background is a{x) = 2d—l for all x, 
to which the addition of a single extra particle causes every site in Z*^ 
to topple infinitely many times. 

We remark that an intermediate behavior is possible, when T„ is infi- 
nite for a finite n, but T„ 7^ Z'^ for all n. An example is the background 
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Figure 2. An exploding sandpile started from n = 5000 
particles in Z^. Background height is 2 except for sites 
in the lattice generated by (1, 10) and (10, 1), which have 
background height 3. Unstable sites are colored black. 

with 3 particles at every site on the x-axis in Z^, and 2 particles at ev- 
ery other site. Adding one particle at the origin produces an infinite 
avalanche of topplings, but each site topples only finitely many times. 
This example shows that exploding is a strictly stronger condition than 
having an infinite avalanche of topplings. 

The papers [6] and [ ] investigated the case of a robust constant 
background of h < 2d — 2 particles at every site. In the regime h < d, 
the diameter of T„ grows like n^^'^; the best known bounds can be found 
in [10, Theorem 4.1]. In the case h = 2d — 2, the set T„ is a cube for 
every n, and an upper bound for the radius is n [6, Theorem 4.1]. No 
proof was found for a better upper bound, even though simulations 
clearly indicated a growth rate proportional to n^/'^. 

In this paper we complete the picture by deriving an upper bound of 
order n^^^ on the diameter of for all h < 2d — 2, and even for some 
backgrounds arbitrarily close to 2(i — 1. We first correct a gap in the 
proof of the outer bound of [10, Theorem 4.1], which we thank Haiyan 
Liu for pointing out to us. Then we use this theorem together with a 
new technique of "background modification" to extend the bounds to 
higher values of h. 

Throughout the paper, we will typically use the symbol a to indi- 
cate a stable background configuration, and rj to indicate an arbitrary 
(possibly unstable) configuration. We use h to denote the constant 
configuration a{x) = h, and we denote a single particle at the origin 
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by do- Write 

Qr = {x : max \xi\ < r} 

for the cube of side length 2r + 1 centered at the origin in l/. Let Ud 
be the volume of the unit ball in W^. 

Our main result, proved in section 3, is the following. 

Theorem 1.1. Fix integers d < h < 2d — 2, and let T^^h 
of sites in Z*^ that topple during the stabilization of h + n6o- Then for 
any e > 0, we have 

for all sufficiently large n, where 

d + e f n 

r = I — 

2d — 1 — h \ujd 

In the case d = h = 2, Theorem 1.1 gives a bound of ^a/w ^ 

l.lSi/n on the radius of the square of sites that topple. Large scale 
simulations by David Wilson indicate that the actual radius is approx- 
imately 0.75^/n. 

It is natural to ask what happens when the background height h 
exceeds 2d — 2. While the background 2d — 1 is explosive, our next 
result shows that there exist robust backgrounds in which an arbitrarily 
high proportion of sites have 2d — 1 particles. For m > 1, let 

A(m) = {x G Z'^ : m / for all 1 < i < d}. 

Thus A(m) is a union of cubes of side length m—1. The following theo- 
rem generalizes the case h = 2d — 2of Theorem 1.1, which corresponds 
to the case m = 1. 

Theorem 1.2. For any m > 1, the background 

a = 2d — 2 + lA{m) 

is robust on Z'^. Moreover, writing Tn^d,cr for the set of sites in Tj^ that 
topple during the stabilization of a + n5o, then for any e > we have 

for all sufficiently large n, where 

(Tl 


On the basis of this theorem, one might guess that 2(i — 1 is the 
critical density below which a background is robust and above which 
it is explosive. Our next two results show that this is not the case. 
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Starting from background height 2d — 2, we can destroy robustness by 
adding extra particles on an arbitrarily sparse lattice L (ZTL^ (Figure 2). 

Proposition 1.3. Let Xj = (xji, . . . , Xj^) jor i = 1, . . . ,d be linearly 
independent vectors in Z'^ satisfying gcd(xij, . . . jX^j) = 1 for all j = 
1, . . . ,d. Let L = Zxi + . . . + ZX(i. Then the background 2d — 2 + 1^ 
on Z"' is explosive. 

Comparing Theorem 1.2 and Proposition 1.3, we see that the geom- 
etry of the extra particles plays a more important role in determining 
robustness than the density of sites at which particles are added. In 
particular, the fact that L intersects every hyperplane in Z'^ that is 
parallel to one of the coordinate hyperplanes, while A(m) does not, 
plays a key role in the proofs. 

As our next result shows, the lattice structure is not essential in 
Proposition 1.3. We can also produce an explosive background by 
adding particles at rare random sites. 

Proposition 1.4. Fixe > 0, and let {l3{x))^(zzd be independent Bernoulli 
random variables with P(/5(a;) = 1) = e. With probability 1, the back- 
ground 2d — 2 + (3 on 7/ is explosive. 

Our proofs make extensive use of the abelian property of the abelian 
sandpile model, which we state and generalize in the next section. 

2. Least Action Principle 

We begin by recalling the notion of toppling procedure defined in [^]. 
This formalism includes most of the natural ways to topple, including: 
discrete time parallel updates, in which all unstable vertices topple 
simultaneously; toppling in nested volumes, in which we successively 
stabilize larger and larger finite regions of Z'^; and Markov toppling in 
continuous time, in which each site has a Poisson clock and attempts to 
topple whenever its clock rings. The technical details of the toppling 
procedures are tangential to our main argument, so the reader may 
wish to skim them and move on to the "least action principle," which 
is the only new material in this section. 

Let X = 77' . We think of elements of X as particle configurations 
on 7/ in which some sites may have a negative number of particles. We 
endow X with the Borel a-algebra coming from the product topology, 
with Z having the discrete topology. On Z"' and on N we use the full 
power set as a cr-algebra, and on the half-line [0, oc) we use the usual 
Borel (T-algebra. A toppling procedure is a measurable function 

T : [0, cx)) X Z"^ X A" ^ N 
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satisfying for all rj E X and all x 

(a) T(0,a;,7;) = 0. 

(b) The function t — > T{t, x, rj) is right-continuous and nondecreas- 
ing with jumps of size at most one, i.e., for all t > 0, 

T{t,x,r]) -limT{s,x,r]) < 1. 

sit 

(c) In every finite time interval, there are only finitely many jumps 
at X. 

(d) There is no "infinite backward chain of topphngs," i.e., no path 
Xi ^ X2 ^ ■ . . and sequence of times ti > t2 > ■ ■ ■ such that for 
all i = 1, 2, . . . 

T{ti,Xi,r]) > \imT{s,Xi,r)). 

We interpret T{t, x, rf) as the number of times x topples in the time 
interval [0, t], for initial configuration 77. We say that x topples at time t 
for initial configuration 77, if 

T{t,x,r)) > \imT{s,x,r)). 

The toppling procedure T may be deterministic (as in parallel updates 
or nested volumes) or random (as in Markov toppling). A random 
toppling procedure can be viewed as a measurable function 

T : [0, 00) X Z'' X A" X Q ^ N U {00} 

where O is a probability space. In this case, we require T(-, -, -,00) to 
satisfy properties (a)-(d) for almost all a; e Jl. 

If u, V arc functions on Z'^, write u < v if u{x) < v{x) for all x £ Z''. 
The discrete Laplacian Au of u is the function 

Au{x) = ^ u{y) - 2d u{x) (1) 

where the sum is taken over the 2d lattice neighbors of x. 

Given a toppling procedure T and initial configuration rj, the result- 
ing configuration at time t is 

rit -.= 7] + AT{t,-,ri). 

We say that T is legal for t] if for all x e Z*^ and alH > such that x 
topples at time t, we have 

lim r}s{x) > 2d. 

sit 

That is, in a legal toppling procedure only unstable sites are toppled. 
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We set 

T(oo, X, rj) = sup T{t, x^rf) G N U {oo}. 

We say that T is /inzte for initial configuration 77, if T{oo,x,ri) < 00 
for all X G Z*^. In this case, we say that T is stabilizing for r/ if 

Voo ■=V + AT(oo, 1]) <2d- 1. 

That is, in the final configuration r/oo no site is unstable. 

The following lemma has been proved a number of times in various 
settings [2, 3, 4, 8, 11]. 

Lemma 2.1. (Abelian property) For any rj & X , if T is a finite legal 
stabilizing toppling procedure for rj, then any legal toppling procedure is 
finite fort]. Moreover ifT' is another legal stabilizing toppling procedure 
for f], then for all x G Z'^ 

T(oo, X, rj) = T'(oo, x, rj). 

If ?7 is a particle configuration for which there exists a finite legal 
stabilizing toppling procedure, then we say that rj stabilizes; otherwise, 
we say that 77 is exploding. If 77 stabilizes, then the function m : Z'' ^ N 
given by 

= T(oo, X, 77), (2) 

where T is any legal stabilizing toppling procedure for rj, is called the 
odometer of rj. 

Note that if every site topples at least once, then rj must be explod- 
ing. Otherwise, by the no infinite backward chain condition (d), some 
site X G Z'' must finish toppling no later than all of its neighbors do; 
since each neighbor topples at least once more, x receives 2d additional 
particles and must topple again, a contradiction. Thus we have shown 

Lemma 2.2. [ , Theorem 2.8, item 4] // rj stabilizes, then u{x) = 
for some x G Z"^. 

Let 7/ be a particle configuration that stabilizes, and let u be its 
odometer function (2). In our application, we will take r] = a + n6o, 
where a is a robust background. Since Au{x) counts the net number 
of particles exiting the site x, the stabilization 7/00 of rj is given by 

Voo=V + ^u. 

Definition. Given a particle configuration rj on l/, a function Ui : 
7/ — > Z, is called stabilizing for rj if 



r/ + Ami < 2rf- 1. 
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Informally, we may think of 77 + Aui as the configuration obtained 
from 7] by performing ui{x) topplings at each site x G Z'^. Note, how- 
ever, that the above definition makes no requirement that these top- 
plings be legal; that is, they may produce sites with a negative number 
of particles. 

Our proof of Theorem 1.1 rests on the following lemma, which char- 
acterizes the odometer function u as minimal among all nonnegative 
stabilizing functions. Deepak Dhar has aptly called this a "least action 
principle," in the sense that the number of topplings in a legal toppling 
sequence is the minimum number required to stabilize the configura- 
tion. In fact, more is true: not only is the total number of topplings 
minimized, but each vertex does the minimum amount of work required 
of it to stabilize the configuration. 

According to the abelian property, if we use a legal toppling proce- 
dure to stabilize 77, then each site x topples exactly u{x) times, regard- 
less of the choice of procedure. The least action principle says that 
in any sequence of topplings that stabilizes 1], even if some of those 
topplings are illegal, each site x topples at least u{x) times. 

Lemma 2.3. (Least Action Principle) Let t] be a particle configuration 
on Z'^ that is not exploding, and let u he its odometer. If Ui : Z*^ — N 
is stabilizing for t], then Ui>u. 

Proof. To compare Ui to the odometer, we use the following discrete 
time legal toppling procedure T' . Enumerate the sites in 7/. Call a 
site X E 'L'^ ready if it has at least 2d particles and has toppled fewer 
than ui{x) times. At each time step, if there are any ready sites, topple 
the smallest ready site. 

Write u'{x) = T'{oo,x,r]) for the number of times x topples during 
this procedure. We will show that u' = u. If rj' = rj + Au' is stable, 
then T' is stabilizing as well as legal, so u' = u by the abelian property. 
Otherwise, rj' has some unstable site y. We must have u'{y) = ui{y); 
otherwise, y would still be ready. Writing u" = ui — u', we obtain 

[r] + Aui){y) = r]'{y) + Au"{y) > r,\y) > 2d 

since u"{y) = 0. This contradicts the assumption that ui is stabilizing. 

□ 

We pause here to record a closely related fact. If Mi,M2 are functions 
on Z'^, write min(Mi,M2) for their pointwise minimum. If x G Z'^ is a 
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site where ui{x) < U2{x), then 



< 'y^^Ui{y) — 2dui{x) = Aui{x). 



y^x 



Likewise, if > U2{x), then A min('Ui, M2)(x) < /S.U2{x). So 

Amin(ui,M2) < max(Aui, Am2)- 
As a consequence, we obtain the following. 

Lemma 2.4. If ui andu2 are stabilizing for a , then mm{ui,U2) is also 
stabilizing for a. 



Proof, cr + A min(Mi, M2) < max{a + Aui, a + A'U2) < 2d — 1. □ 

The set of stabilizing functions is also closed under adding any con- 
stant function, giving it the structure of a module over the tropical 
semiring (Z, min, +). A related module is studied in [9]. 

3. Growth Rates 

Fix an integer h < 2d — 2, and let 77 be the configuration h + n6o 
on Z'^. Let Sn be the set of sites that ever topple or receive a particle 
during the stabilization of rj (in [Id] these were called "visited" sites). 
Note that if y receives a particle, then one of its neighbors must have 
toppled. Thus Sn is related to the set T„ of sites that topple by 

Sn = TnU dTn 

where for A C Z'^ we write 

dA = {y eZ'^ : y ^ A,3z e A,z ^ y}. 

Write |x| = (xf + . . . + x'^Y^'^ for the Euclidean norm on Z*^, and 
for r > let 

Br = {xeZ'^: \x\ < r} 

be the ball of radius r centered at the origin in Z*^. Let Ud be the 
volume of the unit ball in M"^. For the proof of Theorem 1.1, we take 
starting point the following result of [ ]. 

Theorem 3.1. [ , Theorem 4.1] Fix an integer h < 2d — 2. For any 

n > 1, we have 

where r is such that n = UdV'^, and 

ci = (2c/ - 1-/1)- 1/'^ 
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and C2 is a constant depending only on d. Moreover if h < d — 1, then 
for any n> 1 and any e > we have 

Sn C (3) 

where 

c[ = (d-e-/,)-Vrf 
and C2 is independent of n but may depend on d, h and e. 

Note that h may be negative, in which case the background h corre- 
sponds to each site in Z'^ starting with a "hole" of depth H = —h. 

We are grateful to Haiyan Liu for pointing out a gap in the proof 
of the outer bound (3). The gap occurs in Lemma 4.2 of [10], which 
is valid only for H > 0. We correct this gap in section 3.1. Next, 
in section 3.2, we explain our technique of "background modification," 
and use it to deduce Theorem 1.1 from Theorem 3.1. 

3.1. Low Background Height. Fix < /i < c/ — 1, let 77 be the 
configuration h + n6o on Z'^, and consider the odometer function 

Un{x) = number of times x topples during the stabilization of 77. 

The normalization of the odometer function and of the discrete Lapla- 
cian (1) differs by a factor of 2d from the one used in [10]. It is the 
most convenient normalization for the abelian sandpile, since 2d parti- 
cles move in every toppling. 

In [10] it is proved that for every site x G Z'^ with c'j^r — 1 < |x| < c'^^r 
we have 

Un{x) < C 

where c is a constant which may depend on d, h and e but not on n. 
(In the notation of [ ], c = c'2/2d.) 

It follows that Un is uniformly bounded outside the ball Bc'y, indeed, 
if I a; I > c[r, then setting n' = ^ujd{\x\/c[)'^~\, since n < n' we have by 
the abelian property 

Un{x) < Un'{x) < C. 

The next lemma shows that in fact, = outside the slightly larger 
ball Bc'^r+c-i- Hence T„ C Bc'^r+c-i, and hence Sn C i^c^r+o which 
completes the proof of (3). 

Lemma 3.2. For all j = 0, 1, . . . , c and allx G Z'^ with \x\ > c[r+j — l, 
we have 

Un{x) <c-j. 
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Proof. Let 

Rj = {x e Z'^ : \x\ > c[r 

Note that for any x & Rj, all neighbors y ^ x lie in Rj-i, and at least d 
neighbors have \y\ > \x\, so at least d neighbors lie in Rj. 
We will prove the lemma by induction on j. Let 

Uj = {x E Rj I Un{x) > C — j}. 

If Uj is empty, then the proof is complete. Otherwise, by the no infi- 
nite backward chain condition, there exists a site x G Uj that finishes 
toppling no later than all of its neighbors in Uj. By the inductive 
hypothesis, every neighbor y of x satisfies. 

Uniy) <C-j + l. 

Just before x topples for the last time, each neighbor y G Uj has not yet 
toppled for the last time, so y has toppled at most c—j times. Moreover, 
each neighbor y E Rj — Uj has toppled at most c — j times; and each 
neighbor y ^ Rj has toppled at most c — j + 1 times. Hence, just before 
it topples for the last time, x has received at most d{c — j) +d{c — j + 1) 
chips and emitted at least 2d{c — j) chips, leaving it with at most h + d 
chips. Since h < d — 1, this is not enough chips to topple, which gives 
the required contradiction. □ 

3.2. High Background Height. To prove Theorem 1.1 using the 
least action principle (Lemma 2.3), for each coordinate i = 1, . . . ,d we 
will construct a toppling function gi supported in the slab 

Ai^r = {xeZ"^ : \xi\ < r}. (4) 

The effect of toppling according to Qi will be to modify the constant 
background height h by "clearing out" particles down to height at most 
(i — 1 in a smaller slab Ai^^o ^.nd "piling them up" to height at most 
2d — 1 outside Ai,ro- We will see that this can be done while keeping 
To proportional to r. 

On this modified background, n particles at the origin will spread 
with a growth rate at most according to h = d — 1, provided n is small 
enough so that the particles do not spread outside Qr^. This growth 
rate is controlled by Theorem 3.1: n particles on constant background 
height d — 1 in Z'^ spread at most a distance of order n^^^. Since tq is 
proportional to r, we can therefore choose n proportional to r'^. 

The desired background modification can be accomplished a function 
of just one coordinate, gi{xi, . . . ,Xd) = g{xi). The next lemma spares 
us the need to specify g explicitly; it suffices to specify how the back- 
ground is modified. In the lemma, g plays the role of toppling function 
on Z, and / represents the net change in height of the configuration. 
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The conditions (5) mean in words that topphngs cannot change the 
total number of particles, nor the center of mass of a configuration. 

Lemma 3.3. // / : Z Z supported on a finite interval I = [—a,b], 
and 

yel y&I 

then f = Ag for an integer-valued function g supported on the interval 
I' = [1 — a, 6 — 1]. Moreover, if there are no Xi < X2 < such that 
f{xi) < 0, f{x2) > and /(xs) < 0, then g >0. 



Proof. Let 



x-l 



9{x) = {x-y)f{y). 

y=-a 

Then for x > 6 we have 

6 6 

g(x) = X ^ f{y) - Y yf{y) = 

y=-a y=~a 

SO g is supported on I'. Also 

Ag{x) = g{x + 1) - 2g{x) + g{x - 1) 

x-2 



fix) + 5^ ((x + 1 - y) - 2(a; - y) + (x - 1 - y))f{y) 



y=-a 



fix) 



as desired. 

If g{z) < for some z, then since g{—a) = g{b) = 0, the difference 
Dg{y) = g{y + 1) - giy) 

satisfies Dg{yi) < and Dg{y2) > for some yi < z < y2. Hence the 
second difference 

fix) = Ag{x) = Dg{x) - Dg{x - 1) 

satisfies /(xi) < 0, /(X2) > and /(xs) < for some xi < yi < X2 < 
y2 < X3. □ 

Proof of Theorem 1.1. For each i = 1, ... ,c/ we will construct a non- 
negative function Ui on Z'' which is stabilizing for the configuration 
h + n6o, and supported on the infinite slab Aiy, see (4). 
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By the least action principle, Lemma 2.3, the odometer function u 
satisfies u < Ui for i = 1, ... ,d. Since T^^ih the support of u, we 
obtain 

d 

'^n,d,h — f^^i,r = Qr- 

i=l 

To construct Uj, let w : Z*^ ^ N be the odometer function for the 
configuration d — 1 + n6o. By Theorem 3.1, if n is sufficiently large, 
then w is supported on the ball centered at the origin of radius 

In particular, w vanishes outside the cube Qp. 

Let ro be the smallest integer multiple of2d—l — h exceeding p, and 

let 

= 2d-l-h'''- 

Let / : Z ^ Z be given by 

'2{d-l-h), x = 

„. , d—l — h, < |a;| < To 
j{x) = < 

2d — 1 — h, To < \x\ < ri 

^0, |a;| > ri. 

Then with / = [1 — ri, ri — 1] 

J2 fiy) = 2ro(rf -l-h) + (2ri - 2ro)(2ci - 1 - h) 

= -2dro + 2{2d - 1 - h)ri = 0. 

Since f{y) = f{-y) we have Y^yaiVfiy) = 0- By Lemma 3.3, f = Ag 
for a nonnegative integer-valued function g supported on the interval 
J'= [2-ri,ri-2]. 

For X = (xi, . . . , Xd) € Z'^, define 

Ui{x) = w{x) + g{xi). 

Note that the function g{xi) has (i-dimensional Laplacian f{xi). Inside 
the cube Qp, we have /(xj) < d — 1 — h, hence inside Qp 

h + n6o + Aui < d - l + n6o + Aw <2d-l. 

Outside Qp, since w vanishes and /(xj) < 2d — 1 — h, we have 

h + n6o + Aui <2d-l. 



14 



ANNE FEY, LIONEL LEVINE, AND YUVAL PERES 




Figure 3. Left: The stable configuration 2 + n6o + 
Amin(Mi,M2) constructed in the proof of Theorem 1.1. 
Sites with negative height, along the diagonals of the 
square, are colored orange. Right: The stabilization of 
2 + n6o. Here n = 10^. 



Thus Ui is stabilizing for h + n5o- Moreover, since 
ri<d' ^ 



2d-l-h 

l/d 

+ d 



d + e/2 f n 



2d — 1 — h \uOd 

we have ri < r for sufficiently large n, hence Ui is supported on the 
slab Ai^r as desired. □ 

We remark that in addition to bounding the set of sites T^^j^ that 
topple, the proof gives a bound on the odometer function 

u{x) = # times x topples in the stabilization of /i + n6o in Z"^, 

namely 

u{x) < min(Mi(a;), . . . ,Ud{x)) 

= w{x) + mm{g{xi), g{xd)) 
= w{x) + g{max 

By Lemma 2.4, the right side is stabilizing for h + n6o- The resulting 
stable configuration in the case d = h = 2 is pictured in Figure 3. 
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The proof of Theorem 1.2 is identical to that of Theorem 1.1, except 
that we now choose 



/(^) = 



2 - 


2d, 


x = 


—d, 




< X < ro and m \ x 


1 - 


d. 


< X < Tq and m\x 


< 

0, 




fo ^ 2- < ''"i m\ X 


1, 






0' 




\x\ > Ti, 



where ro the smallest integer exceeding P=(l + ^)(;^) , and 

ri = m{drQ — 1). 

Then we have that 

ri < m ^c/ + ^-^^ +m{d-l), 
so that again ri < r for sufficiently large n. 

4. Robust and Explosive Backgrounds 

Write tp^ = ei,...,^d = e^, ipd+i = -ei, ...,^2d = -e^ for the 2d 
coordinate directions in Z'^. If is a rectangular prism in Z^, write 

FiiR) = {y:y^R,y-i;,eR} 

for the outer face of R in direction ipi. 

We will deduce Propositions 1.3 and 1.4 from the following slightly 
more general result. 

Theorem 4.1. Let a be a background on Z'^ satisfying 

(i) a{x)>2d-2 for all x e Z'^; anc? 

(ii) There exists ro G N such that for all r > r^, each outer face 
Fi{Qr) contains a site x with a{x) >2d — 1. 

Then a is explosive. 

Proof. By Lemma 2.2, in order to prove that a configuration rj on 7/ 
is exploding, it suffices to find a toppling procedure in which every site 
in 7/ topples at least once. 

From [(i. Theorem 4.1], if the background height is exactly 2d — 2, 
then for every n, the set of sites that topple during stabilization forms 
a cube Qr, and we can choose n so that r > tq. 

Let -Ro = Qr and 

Rk = Rk-l U -Ffc mod 2d{Rk-l), k > 1. 
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We will define a toppling order in stages /c = 0, 1, 2, . . . so that at the 
end of stage k, all sites in Rk have toppled at least once, and no other 
sites have toppled. Since Ufe>o -^^ ~ follows that every site in Z*^ 

topples at least once, so a + n6o is exploding. 

During stage 0, we perform all the topplings that occur in the sta- 
bilization of 2(i — 2 + n6o- Then every site in the cube Rq has toppled 
at least once, and no other sites have toppled. Hence by (i), every site 
in every outer face of Rq now has at least 2d — 1 particles, and by (ii), 
in every outer face there is at least one unstable site. 

The remaining stages are defined inductively. After stage k — 1, every 
site in F = F/j mod 2d{Rk-i) has at least 2d — 1 particles, and at least 
one site y & F is unstable. Topple first y, then its neighbors in F, then 
the sites in F at distance 2 from y, and so on, until all sites in F have 
toppled once. Now every site in R^ has toppled at least once, and no 
sites outside Rk have toppled, completing the inductive step. □ 

To deduce Proposition 1.3 from Theorem 4.1, since 

gcd{xij,...,Xdj) = 1, j = l,...,d 

there exist integers aij for 1 < i, j < d, satisfying 

d 

^ ^ O-ij^ij 1) j 1, . . . , d. 

Then for each j = 1, . . . ,d, the vector 

d 
i=l 

has Cj-coordinate vjj = 1, so any hyperplane in Z'^ parallel to one 
of the coordinate hyperp lanes intersects L. Moreover, L contains the 
vectors Dej for j = l,...,d, where D = det{xij)fj^i ^ 0. Thus L 
intersects every face Fi{Qr) when r > \D\/2. 

To deduce Proposition 1.4, it remains to check that the configuration 
2d — 2 + f] on Tj'^ satisfies condition (ii) of Theorem 4. 1 with probabil- 
ity 1. Write Si^r for the event that P{x) = for all x G Fi{Qr). By the 
independence of the Bernoulli random variables this event has 

probability 

P(^^_^) = (l_e)l^.(OOI < (1-e)-. 

In particular, Xlr>i Yld=i ^{^i,r) < oo- By the Borel-Cantelli lemma, 
with probability 1 only finitely many of the events £i^r occur. We 
remark on the similarity between this argument and Straley's argument 
for bootstrap percolation [5]. 



GROWTH RATES AND EXPLOSIONS IN SANDPILES 



17 



We define the box Br as 

2d 

Br = dQr = \jF,{Qr). 
i=l 

The following theorem is a partial converse to Theorem 4.1, and gives 
a counterexample to Remark 7.1 in [7]. 

Theorem 4.2. Let ri,r2, ... be an increasing sequence of positive in- 
tegers. Let a be a stable background on in which every site in 
Bri U Br2 U . . . has at most 2d — 2 particles. 
Then a is robust. 

The proof uses the following lemma, which follows from [6, Lemma 4.2] 
and the abelian property. 

Lemma 4.3. If a is a stable background satisfying a{x) <2d — 2 for 
all X & Br, then no sites outside Qr topple during the stabilization of 
a + 5o. 

Proof of Theorem 4-2. We need to show that a+ndo stabilizes in finitely 
many topplings, for every n G N. We induct on n to show that no sites 
outside Qr^ topple during the stabilization of a + ndo- 

By Lemma 4.3, no sites outside Qn topple during the stabilization 
of 0" + (Jo- 
Let (Jn be the stabilization of cr + n5o. By the inductive hypothesis, 
no sites topple outside Qr„ during this stabilization, so an{x) <2d — 2 
for all X G Br^^^. By Lemma 4.3, no sites outside (5r„+i topple during 
the stabilization of o"„ + 5o. By the abelian property, a site topples 
during the stabilization of cr + (n + l)5o if and only if it topples during 
the stabilization of a + n5o or during the stabilization of cr„ + 5o. This 
completes the inductive step. □ 

Remark. Theorem 4.2 remains true for arbitrary disjoint rectangular 
boxes surrounding the origin; they need not be cubical or centered at 
the origin. 

5. Dimensional Reduction 

Our argument used properties of the one-dimensional sandpile to 
bound the growth rate of higher-dimensional sandpiles. There appears 
to be a deeper relationship between sandpiles in d and d—1 dimensions, 
which we formulate in the following dimensional reduction conjecture. 
For X G Z*^, let (Jn,d{x) be the final number of particles present at x in 
the stabilization oi 2d — 2 + ndo- Write rad(n, d) for the radius of the 
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Figure 4. Left: A two-dimensional slice through the 
origin of the sandpile of n = 5 ■ 10^ particles in Tl? 
on background height /i = 4. Right: The sandpile of 
m = 47465 particles in T? on background height h = 2. 
Color scheme on left: sites colored blue have 5 particles, 
turquoise 4, yellow 3, red 2, gray 1, white 0. On right: 
blue 3 particles, turquoise 2, yellow 1, red 0. 

cube {x G Z'^'lan^ix) > 0}. Note that by Theorem 1.1, rad(n, c?) has 
order n^^^. 

Identifying Z^~^ with the hyperplane Xd = in Z'^, we would like to 
compare the slice through the origin Qrad(n,d)nZ'^^^ of the d-dimensional 
sandpile started from n particles with a. {d — l)-dimensional sandpile 
started from some number m of particles. Given m, n and d, let us call 
a site X = (xi, . . . , Xd-i) G Z'^"^ an exact match if 

Crn,d{Xi, Xd-l, 0) = 2 + am,d-l{Xi, Xd-l). 

Given < A < 1, consider the subset Ax of the slice through the origin 

Ax = (<5rad(n,d) " Q\ra.d{n,d)) fl Z"^ ^. 

Conjecture 5.1. There exists a constant X = Xd < 1 such that for all 
n>l there exists m > 1 such that all but 0(rad(n, c?)'^"^) sites in Ax 
are exact matches. 

The case c? = 3 is illustrated in Figure 4. Amazingly, except in a re- 
gion near the origin, the two pictures shown in the figure agree pixel for 
pixel. For some rare values of n, certain small "defects" or "filaments" 
in the two pictures fail to match exactly, which is why we exclude up 
to 0(rad(n, sites. For simplicity, we have restricted our formu- 

lation to the slice through the origin, but dimensional reduction seems 
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to occur in all slices except for those close to the boundary of the cube. 
The value of m is the same for all of these slices. We first learned of 
the dimensional reduction phenomenon from Deepak Dhar. 
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